BCS-BEC crossover at finite temperature in spin-orbit coupled Fermi gases 
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By adopting a T-matrix-based method within the GoG approximation for the pair susceptibility, we study the 
effects of the pairing fluctuation on the three-dimensional spin-orbit coupled Fermi gases at finite temperature. 
The critical temperatures of the superfluid/normal phase transition are determined for three different types of 
spin-orbit coupling (SOC): (1) the extreme oblate (EO) or Rashba SOC, (2) the extreme prolate (EP) or equal 
Rashba-Dresselhaus SOC, and (3) the spherical (S) SOC. For EO- and S-type SOC, the SOC dependence of the 
critical temperature signals a crossover from BCS to BEC state; at strong SOC limit, the critical temperature 
recover those of ideal BEC of rashbons. The pairing fluctuation induces a pseudogap in the fermionic excitation 
spectrum in both superfluid and normal phases. We find that, for EO- and S-type SOC, even at weak coupling, 
sufficiently strong SOC can induce sizable pseudogap. Our research suggests that the spin-orbit coupled Fermi 
gases may open new means to the study of the pseudogap formation in fermionic systems. 

PACS numbers: 03.75.Ss, 05.30.Fk, 67.85.Lm, 74.20.Fg 



I. INTRODUCTION 

The experimental realization of ultracold Fermi gases with 
tunable interatomic interaction has opened new era for the 
study of some longstanding theoretical proposals in many- 
fermion systems. One particular example is the smooth 
crossover from a Bardeen-Cooper-Schrieffer (BCS) super- 
fluid ground state with largely overlapping Cooper pairs to 
a Bose-Einstein condensate (BEC) of tightly bound bosonic 
molecules — a phenomenon suggested many years ago [ 1- 
4]. For a dilute Fermi gas in three dimensions (3D) with a 
short-range interatomic interaction where the effective range 
r of the interaction is much smaller than the interatomic dis- 
tance, such a BCS-BEC crossover can be characterized by 
the dimensionless gas parameter, l/(fcpa), where kp is the 
Fermi momentum and a is the s-wave scattering length of the 
short range interaction. The BCS-BEC crossover occurs when 
l/(&Fa) is tuned from negative to positive values. 

This BCS-BEC crossover has been successfully demon- 
strated in ultracold Fermi gases where the s-wave scatter- 
ing length is tuned by means of the Feshbach resonance [5- 
7]. This has been regarded as one of the key successes in 
the cold atoms researches. Curiously, it is always interest- 
ing to look for other mechanisms of realizing the BCS-BEC 
crossover. Recent experimental breakthrough in generating 
synthetic non-Abelian gauge field in bosonic gas of 87 Rb 
atoms has opened the opportunity to study the spin-orbit cou- 
pling (SOC) effects in cold atomic gases [8]. For fermionic 
case, some of the recent theoretical results suggested that tun- 
ing the SOC may provide an alternative way to realize the 
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BCS-BEC crossover [9-18]. The experimental exploration of 
the spin-orbit coupled Fermi gases has also achieved remark- 
able progresses: The spin-orbit coupled 40 K and 6 Li atoms 
have been realized at Shanxi University [19] and at Mas- 
sachusetts Institute of Technology (MIT) [20], respectively. 

The SOC of fermions can be induced by a synthetic uni- 
form 5/7(2) gauge field, A% — A,<5f , where Aj will play the 
roles of the SOC strengths. With this gauge field, the single- 
particle Hamiltonian reads ft = k 2 /(2m) — er • s(k) where 
s(k) = (X x k x , Xyky, X z k z ). In a very interesting paper [21], 
Vyasanakere and Shenoy studied the two-body problem of 
this Hamiltonian. They paid particular attention to three spe- 
cial types of gauge field configurations: (1) A x = X y = and 
X z = X (called extreme prolate (EP)), (2) X x = X y = X and 
X z = (called extreme oblate (EO)), and (3) X x = X y = X z 
( called spherical (S)). The EO SOC is physically equivalent 
to the Rashba SOC which has been famous in condensed mat- 
ter physics. The EP SOC is physically equivalent to an equal 
mixture of Rashba and Dresselhaus SOCs. The most surpris- 
ing finding of Vyasanakere and Shenoy was that for EO and S 
SOCs, even for a < where the di-fermion bound state can- 
not form in the absence of SOC, the di-fermion bound state 
(referred to as rashbon) always exists and its binding energy 
is generally enhanced with increased SOC. Meanwhile, the 
bound state also possesses non-trivial effective mass which is 
generally larger than twice of the fermion mass m. For two 
dimensional (2D) case, although a bound state exists for arbi- 
trarily small attraction, it was shown in Ref . [ 1 6] that the EO 
or Rashba SOC can generally enhance the binding energy and 
the effective mass of the bound state. 

The novel bound state emerged in the two-body problem 
suggests that the EO or S SOCs may trigger a new type of 
BCS-BEC crossover in the many-body problem of fermions. 
In fact, theoretical studies revealed that for EO or S SOC, even 
at small negative kpa, a crossover from the BCS superfluid to 
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the BEC of rashbons can be achieved by tuning the SOC A 
to large enough value [9-18]. It was shown that for EO or S 
SOCs the system enters the rashbon BEC regime at A ~ vp 
where v-p = k-p /m is the Fermi velocity. Similar conclusions 
were also found for 2D Fermi gases with EO SOC [ 16]. 

So far, most of the theoretical studies of the BCS-BEC 
crossover in 3D SOC Fermi gases focused on the zero- 
temperature ground state based on mean-field theory (MFT). 
Although the MFT captures some qualitative features of the 
zero-temperature crossover, it loses the effects of the pair- 
ing fluctuation which becomes substantial when the system 
goes toward finite-temperature and/or the BEC regimes. In 
the absence of the SOC, previous theoretical studies have al- 
ready revealed that, as a consequence of the pairing fluctua- 
tion, a "pseudogap" emerges in the fermionic excitation spec- 
trum [22-31]. This pseudogap is negligibly small at BCS limit 
but increases as l/(fcpa) is increased and becomes signifi- 
cantly important on the BEC side. Particularly, the pseudogap 
survives above the superfluid critical temperature T c and leads 
to an exotic normal state that is different from the Fermi-liquid 
normal state associated with the MFT. Recently, the experi- 
mental observation of pairing pseudogap in 2D Fermi gases is 
reported [32]. Similar pseudogap phenomena may also appear 
in other strongly correlated systems, such as high-T c super- 
conductors [22, 33, 34], low-density nuclear matter [35], and 
color superconducting quark matter [36]. 

In this paper, we study the spin-orbit coupled Fermi gases 
at finite temperature. To include the pairing-fluctuation ef- 
fects and investigate the possible pseudogap phenomena, we 
will adopt a T-matrix formalism based on a GqG approxi- 
mation for the pair susceptibility which was first introduced 
by the Chicago group [23-28, 33]. This formalism general- 
izes the early works of Kadanoff and Martin [37] and Pat- 
ton [38], and can be considered as a natural extension of the 
BCS theory since they share the same ground state. Moreover, 
this formalism allows quasi-analytic calculations and gives a 
simple physical interpretation of the pseudogap emergence. It 
clearly shows that the pseudogap is due to the incoherent pair- 
ing fluctuation. Within this formalism, we can also determine 
the superfluid critical temperature and study how the pairing 
fluctuation affects the thermodynamics. We note that this is 
the first systematic study of the 3D spin-orbit coupled Fermi 
gases at finite temperature. For 2D spin-orbit coupled Fermi 
gases, the possible BKT transition at finite temperature was 
already studied [16]. 

The article is organized as follows. In Sec. II, we present a 
detailed theoretical scheme of the T-matrix based formalism 
at finite temperature. The numerical results are given in Sec. 
III. We summarize in Sec. IV. Throughout this article, we use 
natural units h = ks = 1. 



II. T-MATRIX-BASED FORMALISM 

We consider a homogenous Fermi gas interacting via a 
short-range attractive interaction in spin-singlet channel. In 
the dilute limit kpro <C 1 and mXro <C 1 where tq is the ef- 
fective interaction range, this system can be described by the 



following Hamiltonian, 



H = 



+ U / d 3 r ^(r)V>J(r)^(r)V t (r), (2.1) 



where Hq = — V 2 /(2m) — /i is the free single-particle 
Hamiltonian with /i being the chemical potential, H so = 
—i Yli=i ^i a idi is the SOC term, and U < denotes the at- 
tractive s-wave interaction. 

Introduce the four-dimensional Nambu-Gorkov spinor ^I" = 
(ipf, ip±, ipp 4>\) T ■ The (imaginary-time) Green's function of 
the Nambu-Gorkov spinor is given by 

S(t,t) = -(T r *(r,r)*t( ,o)) 
~G{t,t) F{t,v)~ 



J"(r,r) 0(r,r) 



(2.2) 



where T r is the (imaginary-) time-ordering operator and r 
it. It is convenient to work in frequency-momentum space, 



S{K) = 



G{K) F{K) 

Hk) q{k) 



(2.3) 



where K = (fco = iw„,k) with oj n = (2n + 1)ttT (n inte- 
ger) being the Matsubara frequency for fermion. The Green's 
functions have the following properties: 

S(iw„,k) = -[g{-iu n ,-k)] T , (2.4) 
J-(iw„,k) = +[7"(-ia; n ,k)]t, (2.5) 

^(iWn.k) = -[T(-lLU n ,-k)} T , (2.6) 

■F(iw„,k) = -[^(-^„,-k)] T , (2.7) 
0(iw„,k) = +[a(-iw„,k)] t , (2.8) 
S(iw„,k) = +[£Hw„,k)] t . (2.9) 

In the rest of this section, we will introduce the basic 
method of the T-matrix. Our strategy will closely follow 
Refs. [23-28, 33, 35]. The T-matrix we will adopt is defined 
as an infinite series of ladder-diagrams in the particle-particle 
channel by constructing the ladder by one free particle prop- 
agator and one full particle propagator. The T-matrix thus 
enters the particle self-energy in place of the bare interaction 
vertex. The equation that defines the T-matrix, the self-energy 
equation (or gap equation) as well as the number density equa- 
tion form a closed set of equations, and should be solved con- 
sistently. One can view this approach as the simplest gener- 
alization of the BCS theory, which can also be cast formally 
into a T-matrix formalism. Let us start with the BCS theory. 



A. BCS Theory 

The BCS theory is based on the mean-field approximation 
to the anomalous self-energy. We start with the mean-field 
inverse fermion propagator, 



ia y A m { 



-ia y A mf -[g- L (-K)Y 



(2.10) 
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where the anomalous self-energy A m f is chosen as constant 
and can be used as order parameter for superfiuid phase tran- 
sition. The inverse free fermion propagator reads, 



Q Q 1 (K) = iuj n - £ k - £ so (k), 



(2.11) 



with £ k = k 2 /(2m) - /i and £ so (k) = £\ =1 X^k, (X, 
real). By direct doing the matrix inversion, one obtains 



S mf {K) 



FmiiK) G^{K) 



(2.12) 



A ll {K) + ^-B 11 {K), 
A 22 (K) + ^B 22 (K), 

A l2 (K) + ^-B 12 (K) 
A 21 (K) + ^-B 21 (K) 



Its elements are 

Gmf(K) = 
Fmf(K) = 



where we introduced A|k| = yY^i=i ^i^i anc ^ 

A n (K) = 

A 22 {K) = 

Ai 2 {K) = 
A 21 {K) = A 12 (K) 



(2.13) 
(2.14) 
,(2.15) 
(2.16) 



1 




iu n + £ k 


2 




^ n ) 2 -(i?k) 2 J 


1 




iu n ~ £ k 


2 




^„) 2 -(£ k ) 2 J 


1 


Amf 




2 


[(iu^ - (E+y 


^„) 2 -( J B-) 2 J 



(2.17) 



and 



Bu(K) 
B 22 (K) 



iu>„ + £ k 



^ n y-{E+f ^„) 2 -(i?k) 2 



B 12 {K) = - 



iuj n - £ 



( lWn ) 2 - (£+) 2 (iw n )a - (£ k ) 2 



( lW „) 2 - (£+) 2 ^„) 2 - (£") 2 



B 21 (K) = -B 12 (K). 



(2.18) 



Here £± = y^) 2 + A 2 m{ with £ = & ± A|k| is the 
fermion dispersion relation. One can verify that Eqs. (2.13)- 
(2.16) satisfy Eqs. (2.4)-(2.9). 

Then from the standard Green function method, the coupled 



gap and density equations are expressed as 

U 

K 

1 - 2n F (E£) 



2f3V 
UA ml 



2V 



K 

EE 



a=± k 



2EP 



(2.19) 



h E E V(*£) + «) 2 "f(-^)] , 

a=± k 

(2.20) 



where n F (x) = l/[exp (/3x) + 1] is the Fermi-Dirac function, 
e" ,LJ " with — >• + is a convergence factor for the Matsubara 
summation, and fug) 2 = ±(1 + and (w£) 2 = i(l - 

^k/^k) are tne Bogoliubov coefficients. 

The difference between Q^^K) and C/~ f (K) defines the 
mean-field self-energy, 

= -A 2 nf za^o(^)^y (2.21) 
If we define a T matrix in the following form, 

t ml (Q) = -A 2 lf <5(Q), (2.22) 

where Q = (qo = icu v ,q) with w„ = 2vkT (v integer) be- 
ing the boson Matsubara frequency and 5(Q) — fi8 v $8^ (q), 
Eq. (2.21) can be rewritten in a manner of 



i(TyGo(K - Q)ia y .(2.23) 



This shows that in the BCS theory, the fermion-fermion pairs 
contributes to the fermion self-energy only through their con- 
densate at zero momentum, and these condensed pairs are as- 
sociated with a T-matrix or propagator (2.22). 

Furthermore, if we define the mean-field pair susceptibility 

as 



Xmf(O) 



1 



2(3V 



Tr E [Gmf(K)ia y G (K - Q)ia y 

K 



(2.24) 

we can rewrite the gap equation in the superfiuid phase as 

1 + Uxmf (0) = 0, T<T C . (2.25) 

This suggests that if one considers the uncondensed pair prop- 
agator or T matrix to be of the form 



^-pair 



u 



U X mi(Q) 



(2.26) 



then the gap equation is given by t~ ir (0) = 0. 

It is well known that the critical temperature T c in the BCS 
theory is related to the appearance of a singularity in a T ma- 
trix in the form of Eq. (2.26) but with A m f = 0. This is the 



4 



so-called Thouless criterion for T c [39]. But the meaning of 
Eq. (2.25) is more general as stressed by Kadanoff and Mar- 
tin [37]. It states that under an asymmetric choice of XmU me 
gap equation is equivalent to the requirement that the T ma- 
trix associated with the uncondensed pairs remains singular at 
zero momentum and frequency for all temperatures below T c . 

Although the construction of the uncondensed pair propa- 
gator (2.26) in BCS scheme is quite natural, the uncondensed 
pair has no feedback to the fermion self-energy (2.23). In the 
BCS limit (both \ U\ and A are small), such a feedback may not 
be important, but if the system is strongly coupled (for large 
\U\ and/or for large A for EO or S SOC), this feedback will be 
essential. The simplest way to include the feedback effects is 
to replace t m f in Eq. (2.23) by t m i -Mpair. But to make such an 
inclusion self-consistent, i pa i r should be somewhat modified 
which we now discuss. 



pair is given by t~^(Q) — 0, hence i pg (0) oc |U& with //{, the 
effective chemical potential of the pairs. Then the BEC con- 
dition requires ^ = 0, and thus t~^(Q) = 0, for all T < T c . 

The gap equation (2.30) tells us that t pg (Q) is highly 
peaked around Q = 0, so we can approximate £ pg as 



£ pg (iT) ~ -A 2 pg ia y g (K)ia y , T<T C , (2.31) 
where we have defined the pseudogap parameter A pg via 



(2.32) 



The total self-energy now is written in a BCS-type form 



B. G G Formalism at T < T c 

The BCS theory involves the contribution to the self-energy 
from the condensed pairs only, but generally, in superfluid 
phase, the self-energy consists of two distinctive contribu- 
tions, one from the superfluid condensate, and the other from 
thermal or quantum pair fluctuations. Correspondingly, it is 
natural to decompose the self-energy into two additive terms, 



= E rf (K) + EpgCA-), 
with the T matrix accordingly given by 



(2.27) 



KQ) 

t si {Q) 
*pg(Q) 



t s{ (Q) + t pg (Q), 

-A S V(Q), 
U 



i + UxiQY 



(2.28) 



where the subscript sf and pg indicate that these terms are 
responsible to the superfluid condensate and pseudogap in 
fermionic dispersion relation. See Fig. 1 for the Feynman 
diagrams for t pg (Q) and 'E(K). Comparing with the BCS 
scheme, t m f(Q) in Eq. (2.23) is replaced by t{Q), and T,(K) 
now contains the feedback of uncondensed pairs. Inspired by 
Eq. (2.24), we now choose the pair susceptibility x(Q) to be 
the following asymmetric GqG form, 



x(Q) 



i 



2f3V 



J2 0(K)i(r y g (K - Q)ia v . (2.29) 



K 



In spirit of Kadanoff and Martin, we now propose the su- 
perfluid instability condition or gap equation as [extension of 
Eq. (2.25)] 



1 + C7 X (0) = 0, T<T C 



(2.30) 



We stress here that this condition has quite clear physical 
meaning in BEC regime. The dispersion relation of the bound 



E(JsT) = -A 2 ia y g (K)ia 



(2.33) 



- Ap g . It is clear that A pg also contributes 



but with A 2 = A 2 

to the energy gap in fermionic excitation. Physically, the 
pseudogap A pg below T c can be interpreted as extra contribu- 
tion to the excitation gap of fermion: an additional energy is 
needed to overcome the residual binding between fermions in 
a thermal excited pair to produce fermion-like quasi-particles. 
One should note that A pg is associated with the thermal fluc- 
tuation of the pairs A pg (T) - (A 2 (T)) - (A(T)) 2 [25, 33] 
hence it does not lead to superfluid (symmetry breaking). Be- 
sides, at T = 0, the GqG formalism recovers the BCS theory, 
hence the GqG formalism does not involve quantum fluctua- 
tion. 



+ 



+ 



J nif 



FIG. 1 : Feynman diagrams for the T matrix of non-condensed pairs 
and the fermion self-energy in the GoG formalism. 



With the self-energy (2.33), it is easy to see that the 
gap equation and number equation remain the forms of 
Eqs. (2.19)-(2.20) except the replacement of A m f — > A: 



2v 2-" 



a=± k 



2E£ 



(2.34) 



1 Tl J2e^g(K) 



/3V 



K 



a=± k 



(2.35) 
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The pair susceptibility (2.29) can be calculated as, 

M»E£)-M-C-k) 



q,7,s— ± k 



5o 



q-k 



with go = iw„ and 



7- kq = E 



Affejfa - kj) 
i A 2 |k||q-k| ' 



(2.36) 



(2.37) 



Furthermore, the gap equation (2.30) suggests that we can 
make the following Taylor expansion for x(Q)> 



where Z is a pair wave-function renormalization factor and 
mu is the effective "boson" mass parameter in the i direction. 
A straightforward calculation leads to 



CK,S — ± k k 

<9x(Q) I 



(2.39) 



Z = 



ov X/ X/ 



2 ^ 

a,s— ± k 



(2.40) 



X(Q) = X(0) + 2 ft 



3 1 

Y — 



i=l 



2m 6i 



(2.38) 



Z 



1 d 2 x (Q) 

2 9g? 



1 



Q=0 

E E 



2s 



a,s— ± k 



2E£ A 2 k 2 



A 2 k 2 



a 



Eg \m 



A|k| 

1 
"A|k| 



a , ,. , — a 



A 3 |k| 3 



A? 



jse y A 2 k 2 



> 2 fc[\ n F (^ k ) - »f(-C 



A 2 k 2 J 



t-a 



*3E 



(2.41) 



According to this Taylor expansion, we apply the pole approx- 
imation to the pair propagator or T matrix t pg (Q), 



t PS (Q) 



9o ~ Ei=i ??/( 2m w) 



(2.42) 



We stress here that in general, the small Q expansion of 
tZg(Q) should contain a term oc g 2 . Without this term, 
Eq. (2.42) does not respect the particle-hole symmetry and 
thus can work, in principle, only when the system becomes 
bosonic. At the BCS limit, the system possesses a sharp Fermi 
surface, the pair propagator should asymptotically recover the 
particle-hole symmetry, i.e., the g 2 term should be kept. How- 
ever, since at BCS limit the pseudogap is expected to be very 
small, applying Eq. (2.42) does not bring much quantitative 
difference. Therefore, we will apply Eq. (2.42) to the whole 
crossover region. 

Substituting the number density equation, the parameter Z 
can be expressed as 



Z = 



1 

A 2 



n 

2 2V 



(2.43) 



a=± k 

The expression in the square bracket of the right-hand side is 



nothing but the density of the pairs rib, we thus have rib 
ZA 2 . 

Substituting Eq. (2.42) into Eq. (2.32) leads to 



1 3 

Z yE nB E^/( 2m fc)] 

q i= 



i=l 



Z 11 V 2tt \2 

i=l v 



(2.44) 



where Ub(x) = l/[cxp ((3x) — 1] is the Bose-Einstein func- 
tion and a vacuum term was regularized out. It should be 
stressed that at zero temperature A 2 g = 0, hence the GqG 
scheme yields the BCS ground state. It is also worth noting 



that A 2 g 



-.uncondcnscd 



/Z, and hence A^ f 



-.condensed 



/z. 



Now, Eq. (2.34), Eq. (2.35), as well as Eq. (2.44) are cou- 
pled to determine the total excitation gap A, the pseudogap 
A pg and the chemical potential \x at a given temperature be- 
low T c , and T c itself is determined by the vanishing of A s f . 
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C. G G formalism at T > T c 

Above T c , Eq. (2.30) does not apply, hence Eq. (2.31) no 
longer holds. To proceed, we extend our more precise T <T C 
equations to T > T c in a simplest fashion. We will con- 
tinue to use Eq. (2.33) to parameterize the self-energy but with 
A = A pg , and ignore the finite lifetime effect associated with 
normal state pairs. In the absence of the SOC, it was shown 
that this is a good approximation when temperature is not very 
much higher than T c [24, 26]. The T matrix t pg (Q) at small 
Q can be approximated now as 



t Pg (Q) 



q - O q 



(2.45) 



where f2 q = ^i=i Qi /{^ m bi) — /-**>■ Since there is no conden- 
sation in normal state, the effective pair chemical potential [i b 
is no longer zero, instead, it should be calculated from 

Z^ ee f- 1 (0) = i + x (0) 

= h-w (2 - 46) 

This is used as the modified gap equation. Similarly, above T c 
the pseudogap A pg is determined by 



Z J-l V 2tt "S 



(2.47) 



where Li„(z) is the polylogarithm function. Then Eq. (2.46), 
Eq. (2.47) and the number equation which remains unchanged 
determine A pg , fi and ^ at T > T c . 



III. RESULTS AND DISCUSSIONS 



in the latter case. In the molecular BEC limit, we expect \i < 
and |/i| 3> T c . For temperature around or below T c , we can 
approximate ££/T,E£/T rj oo, and all the equations be- 
come temperature independent. In this limit, the gap equation 
determines the chemical potential while the number equation 
determines the gap. By further expanding the equations in 
powers of A/|/i| and keeping several leading terms, we ob- 
tain some analytical results for various quantities at T < T c . 
(Some of them are already reported in Refs. [9-12, 17, 21].) 
(I) S case. The chemical potential is well given by 



Eb 
"~2~' 



(3.2) 



where Eb is the binding energy determined by the two-body 
problem [9, 17, 21], 



Er 



raX 2 + 



- 1 

4m \ a 



1 



4to 2 A 2 



(3.3) 



The effective pair mass nib coincides with the molecular ef- 
fective mass determined at the two-body level. We have 



2m _ 7 4 ( E B - m\ 2 \ 3/2 2mA 2 
m b 3 3 V E B ) E B 

Other quantities such as A and Z can be evaluated as 

32e F ^ fE B -m\ 2 \ 3/2 



A 1 



3tt E b V 
n Er 



3/2 



2 A 2 8tt \E b -mX 2 ^ 
(II) EO case. The chemical potential is also given by 



Eb 



(3.4) 

(3.5) 
(3.6) 

(3.7) 



where two-body binding energy Eb is determined by the al- 
gebra equation [9-11,21], 



With all the equations settled down, we now present the pre- 
dictions obtained by solving them numerically. We will focus 
on three different types of SOC: EP (A^ = X y = 0, X z = A), 



EO (A, 



A,, 



A, X z = 0), and S (X x = X y = X z = A). 



In all these cases, we regularize the UV divergence in the 
gap equations by introducing the s-wave scattering length a 
through 



1 

U 



m 
4:T:a 



v 4^ 



m 
k 2 ' 



(3.1) 



A. Analytical Results in the Molecular BEC limit 

Let us first examine the molecular BEC limit which can be 
achieved by either tuning l/(fcpo) —> +oo for fixed SOC or 
tuning A — s- oo for fixed l/(fcpa) for EO and S SOCs. The 
former case is well studied and here we are mainly interested 



E B _ 1 V^+x/^A 2 
mA 2 2 n JEb - Vm~X^ 



T-- (3.8) 
Xa 



The effective pair mass becomes anisotropic and is given by 



,11 



2m 
2m. 



1 - 



,A 2 



E B — mA 2 , Eb — mX 2 
■ In 



2Eb 2Eb Eb 

% ~ ^"<- (3-9) 
Other quantities such as A and Z can be evaluated as 

8£ F (E B -mX 2 ) [2e^ 



A 2 

Z 



2A 2 



3tt V E B 

\fn¥E~B 
= 8tt(E b - mA 2 )' 



(3.10) 
(3.11) 



(III) EP case. We find that the EP case is trivial. Increasing A 
can not induce a BCS-BEC crossover. For large and positive 
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1/ikFo), the EP SOC only induce a shift for the chemical 
potential, 



1 



mA 2 



2ma 2 2 

The pair effective mass is almost isotropic and is given by 



(3.12) 



II 



2m. 



(3.13) 



Other quantities such as A and Z just recover the usual results 
without SOC, 



A' 



8V2e F k F 
37r ma 



Z 



8\/2V 



(3.14) 
(3.15) 



The critical temperature in the BEC limit, Tbec> is deter- 
mined by the number equation 



n B 



1 



V 



k exp 



(3.16) 



where e B (k) = £\ =1 kf/(2m bi ). This leads to T BE c = 
27r[nB/(Vriim6iC(3/2))] 3 ' 2 in three dimensions. Setting 
iib = n/2, we obtain 



Tbec « 0.218e F f 



2m 

m 6 i 



1/3 



(3.17) 



Therefore, in the molecular BEC limit, Xbec is only a func- 
tion of the combined dimensionless parameter rj = l/(mXa). 
For S SOC we have 

6m, i] — > — oo 
[ 2m, i] -> +oo (3.18) 
2.32m, 77-^0 




0.0726e F , V -> -00 

0.218e F , r/^+oo (3.19) 

0.188e F , V -> 0- 



4m, 77 — > —00 

2m, i] -> +00 (3.20) 

2.40m, 77 -> 



0.137e F , -> -00 
7bec = { 0.218e F , tj -> +00 (3.21) 
0.193e F , -> 0. 

As we will see, the above obtained Tbec coincide well with 
our numerical results in Sec. Ill B. 
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FIG. 2: (Color online) The critical temperature T c scaled by the 
Fermi energy epasa function of the gas parameter 1 / (foo) for fixed 
SOC A = and A = vf- Also shown is the pair dissociation temper- 
ature T*. 



B. Superfluid Critical Temperature 

By numerically solving the set of coupled gap, number den- 
sity, and pseudogap equations, we can obtain the superfluid 
order parameter A s f, the pseudogap A pg , and the fermion 
chemical potential \i. The superfluid critical temperature T c is 
determined by the vanishing of the superfluid order parameter 
A s f. The numerical results for T c /ef as a function of the gas 
parameter l/(fc F a) is shown in Fig. 2 for A = and A = vp- 
T c /ev as a function of the SOC A is shown in Fig. 3 for fixed 
l/(fc F a) = —2 and l/(fc F a) = 00. Also shown is the criti- 
cal temperature predicted by the BCS theory, T*/e F , which is 
determined by the vanishing of A m f . It is monotonously in- 
creasing as l/(/c F a) or X/vv increased. The BCS theory loses 
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the pairing-fluctuation effect and does not give reliable critical 
temperature particularly at large 1 / (kpa) or A/ vp where T c is 
mainly determined by the bosonic degrees of freedom. 

For all three types of SOC, we find that T c is a smooth func- 
tion of \/(kpa) and A/up, and the superfluid phase transi- 
tion is always of second order for the whole crossover region 
(see next subsection). Also, it can be seen that T c is not a 
monotonous function of \/{kpa) when A is small: There is a 
local maximum in T c curve around the unitary point. Simi- 
lar local maximum also appears when one uses the Nozieres- 
Schmitt-Rink approach to determine T c in the absence of 
SOC. It may be understood by noticing that the BEC critical 
temperature is increased when repulsive interactions between 
bosons are turned on [40]. We note that for a Rashba spin- 
orbit coupled Fermi gas, the superfluid transition temperature 
has been roughly estimated by approximating the system as 
a non-interacting mixture of fermions and rashbons and has 
been found to increase monotonously across the BCS-BEC 
crossover [11]. 

From the top panels of Fig. 2 and Fig. 3 we see that the 
EP SOC does not affect T c and T*. This is consistent with 
the observation that EP SOC solely does not lead to new 
novel bound state and the fermion excitation gap does not 
change [21]. This can be understood by noticing that the EP 
SOC in Hamiltonian (2.1) can be gauged away by using the 



gauge transformation 



-iraXz 



-0-)- and ipi — y e m ip±, 



resulting only a constant shift in the chemical potential, /i — » 
(i + mA 2 /2. 

For the EO and S SOCs, we observe from Fig. 3 that T c 
is not sensitive to A for A <C vf and A ^> vf, but becomes 
sensitive to A for A ~ vp'- One can identify a BCS-BEC 
crossover solely induced by A around A«i^. This coincides 
with previous zero-temperature studies [9-18]. Although at 
large negative l/(fcpa), T c increases (almost) monotonously 
as A grows, near the resonance, we find that T c is a decreasing 
function of A, in contrast Ref. [41] where the authors predicted 
an increasing T c along with A. We note that for large enough 
A, our result converges correctly to a universal molecular limit 
Tbec (see Sec. Ill A). 

It should be stressed that T* sets a lower bound for the pair 
dissociation temperature above which the pairs essentially dis- 
sociate due to thermal excitations. Previous study in the ab- 
sence of the SOC shows that it is a good approximation to set 
T* on the BCS side or near unitary as the pair dissociation 
temperature [24, 26]. So we refer to T* as the pair dissoci- 
ation temperature in Fig. 2 and Fig. 3. The region between 
T c and T* is a pseudogap dominated window in which a nor- 
mal state is no longer described by the Landau Fermi liquid 
theory. Fig. 2 and Fig. 3 show that even for large negative 
l/(fcpa) the pseudogap dominated region can be sizable once 
the SOC is large. Thus, the spin-orbit coupled Fermi gas may 
provide a new platform to study the formation of pseudogap 
in fermionic systems. 
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FIG. 3: (Color online) The critical temperature T c and the disso- 
ciation temperature T* scaled by the Fermi energy ef as a func- 
tion of the SOC A/uf for fixed gas parameters l/(kpa) — —2 and 
l/(k F a) = 0. 



C. Pseudogap 

In this subsection, we focus on S and EO SOCs because EP 
SOC does not bring qualitatively new features to the temper- 
ature dependence of the pseudogap than the A = case. In 
Fig. 4-Fig. 7, we plot A, A s f, and A pg as well as Abcs ( m 
units of ef, the same below) as functions of temperature. 

The common feature for all these figures are: A s f monoton- 
ically decrease to zero at T c . Below T c , A pg (T) is a mono- 
tonically increasing function from zero at T — where it 
vanishes according to A pg oc T 3 / 4 (see Eq. (2.44)). Above 
T c , Ap g (r) is a monotonically decreasing function from its 
maximum value located at T c . This kind of temperature de- 
pendence clearly shows that pseudogap is due to the thermally 
excited pairs: Below T c when T goes higher more pairs are 
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FIG. 4: (Color online) The temperature dependence of A, A a f , and 
A pg scaled by the Fermi energy for a Fermi gas without SOC at 
fed = —0.5. 
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excited from the condensate and at T c all condensed pairs are 
thermally excited; after that the thermal motion of the pair par- 
ticipators begins to dissociate the pairs and hence A pg (more 
precisely, ZA^ g ) begins to decrease. Although the physical 
pictures are clear, at temperature much higher than T c our for- 
malism may fail since the finite life-time of the pairs, which 
is not included in our formalism, may become important. 

By comparing Fig. 5-Fig. 6 to Fig. 4 and by comparing two 
bottom panels of Fig. 7 to the panel on the top, one can see 
that although the SOC does not modify the general tendency 
of the temperature dependence of the gaps, large SOC signif- 
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FIG. 6: (Color online) The temperature dependence of A, A a f , A pg , 
and Abcs scaled by the Fermi energy at k^a — —0.5 for EO SOC. 



icantly enlarges the pseudogap window in the normal phase. 
Such pseudogap window may be detected by RF spectroscopy 
measurements, which we now turn to study. 



D. RF Spectroscopy 

The radio-frequency (RF) spectroscopy has been proven to 
be very successful in probing the fermionic pairing, quasi- 
particle excitation spectrum, and superfluidity. For an atomic 
Fermi gas with two hyperfine states, | and | \), the RF 
laser drives transitions between one of the hyperfine states 
(i.e., | I)) and an empty hyperfine state |3) which lies above 
it by an energy oj^ (which to be set to zero because it can be 
absorbed into the chemical potential) due to the magnetic field 
splitting in bare atomic hyperfine levels. The Hamiltonian for 
RF-coupling may be written as, 



Hrf — Vb 



d 3 i 



V4(r)^(r)+H.C. 



(3.22) 



where Tp\ (r) is the field operator which creates an atom at the 
position r and Vb is the strength of the RF drive and is related 
to a Rabi frequency lor by Vb = 

Let us now assume that there is no interaction between the 
third state and the spin-up or spin-down states, i.e., there is 
no final state effect. This approximation sounds valid for 40 K 
atoms, where the s-wave scattering length between the spin- 
down state and the third hyperfine state is small (i.e., ~ 200 
Bohr radii) [42]. Within this approximation and taking into 
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FIG. 7: (Color online) The temperature dependence of A, A s f , A pg , 
and Abcs scaled by the Fermi energy at unitary point for A = 
Fermi gas and for EO- and S- SOC Fermi gases. 
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for SO coupled unitary Fermi gases at resonance. 



account that the third state is not occupied initially, the transfer 
strength (integrated RF spectrum) per spin-down atom can be 
written by (Vq — 1), 

r(w) = T7^ E ^(k, £k - w)n F (£ k - w), (3.23) 

where nj_ is the number density of the spin-down fermion and 
An = — (l/ 71 ")! 111 ^^ i s tne spectral function of the spin- 
down state. Note that 



du>r(ui) = 1, 



(3.24) 



because 



"4 



1 f 00 

- / du^2A u 0L,u)n F (u). (3.25) 

* J — OO I. 



In Fig. 8, we present the integrated RF spectra for three 
different SOCs at resonance and at T c , T*, and 2T* 1 . The RF 
spectra are calculated in an idealized manner (see Appendix 
A), i.e., we neglect the final state interaction and width of the 
uncondensed pairs. Taking into account the effects of finite 
width may change the shape of the RF spectra, however, most 
of the qualitative features shall retain. 



1 Theoretical predictions for RF spectroscopy of single spin-orbit coupled 
bound fermion pair and of noninteracting spin-orbit coupled Fermi gas 
have been reported in Ref. [43]. Recently, the RF spectroscopy of equal 
Rashba and Dresselhaus spin-orbit coupled Fermi gases has been studied 
experimentally and theoretically in Ref. [44]. The theory part is based on a 
formalism similar with what we used here. 
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It is seen from Fig. 8 that the RF spectra consist of two con- 
tinuum branches, one positive and another negative. The pos- 
itive branches correspond to the "binding" fermion pairs con- 
tribution, while the negative branches can be regarded as the 
response of the thermal excited quasi-particles with a pseu- 
dogap A pg reflected in the positions of the negative branch 
peaks. With increasing temperature, more quasi-particles are 
excited, leading to a much more pronounced response appear- 
ing in the negative branches. Such a temperature-sensitive 
feature of the RF spectroscopy may provide a useful way to 
experimentally measure the critical temperature and to detect 
the existence of pseudogap in the normal phase. 



IV. SUMMARY 

We have theoretically investigated thermal effects on the 
BCS-BEC crossover of spin-orbit coupled Fermi gases. For 
this purpose, we have employed a T-matrix formalism based 
on a GoG approximation for the pair susceptibility, which was 
thoroughly used in the previous studies of Fermi gases with- 
out SOC. This formalism extends the standard BCS theory by 
appropriately decomposing the excitation gap to a condensa- 
tion part and a pseudogap part that characterize the pairing 
fluctuations. 

Comparing to the BCS theory, our GoG formalism pre- 
dicts lower and more reliable critical temperature T c for the 
superfluid-normal phase transition. The results for T c have 
been presented in Fig. 2 to Fig. 3. At various molecular BEC 
limits, our predictions correctly recover the BEC temperature 
Tbec °f f ree Bose gases. The pseudogap persists not only 



J 



in the superfluid phase (T < T c ) but also in a window of 
the normal phase (T c < T < T*) where it represents the 
existence of non-condensed, preformed pairs which dissoci- 
ate above T* . We have studied how the SOC influences the 
emergence of the pseudogap, as shown in Fig. 4-Fig. 7. It is 
seen that strong S- or EO- type SOC can significantly enlarge 
the pseudogap window in the normal phase. Thus, spin-orbit 
coupled Fermi gases provide a new platform to study the pseu- 
dogap physics. Experimentally, such pseudogap might be re- 
vealed by RF spectroscopy measurements. We have presented 
our qualitative predictions on the RF spectra in Fig. 8, which 
may be easily tested in future experiments. 
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Appendix A: Expressions for the RF spectroscopy 

In this appendix, we list some expressions for the ideal- 
ized RF spectroscopy for Fermi gases with and without SOC. 
Some of them are used in Sec. HID. By "idealized", we mean 
that these expressions neglect the effects due to final state in- 
teractions and due to the finite lifetime effects of the uncon- 
densed pairs. 



(I) Without SOC. The spectral function of spin-down fermion is given by 

A° u (k, u) = ul6(u - E v ) + vlS{cj + E k ) 

The RF spectrum: 

m 3/2 A 2 Ilo 2 - A 2 + 2[iuj 



(Al) 



r (w) = 



47T 2 LO 2 



-np 



2w 



e 



A 2 + 2/mj 



(A2) 



(II) S SOC. The spectral function spin-down fermion is given by 



Ok 11/ 

The integrated RF spectrum: 

1 f°° 

TsH - i^E/ o d l k H k l 8 [(<) 2 S(&-"-EZ) + Kf6(&-oj + EZ)}n F (&-u>). 



(A3) 



(A4) 



(III) EO SOC. The spectral function of spin-down fermion is given by 

A™(k,u) = ^[(v£)H(u,-EZ) + (vZ)*6(u, + ££)]. (A5) 

a 
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The integrated RF spectrum: 
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- A 2 - A 2 fc 2 fc 2 



2(u) + a\k^ 



x9 



w 2 - A 2 - X 2 k 2 , 



2(u + aXk±) 2m 
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(A6) 



(IV) EP SOC. The spectral function of spin-down fermion is given by: 



The integrated RF spectrum: 



r E p(w) 



87r 2 nj, 



/>oo 
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